A CHANDRA STUDY OF TEMPERATURE 
SUBSTRUCTURES IN INTERMEDIATE-REDSHIFT GALAXY 

CLUSTERS 



Liyi Gu^, Haiguang Xu^, Junhua Gu^, Yu Wang^, Zhongli Zhang^, Jingying Wang^, 
Zhenzhen Qin^, Haijuan Cui^, Xiang-Ping Wu^ 

ABSTRACT 

By analyzing the gas temperature maps created from the Chandra archive 
data, we reveal the prevailing existence of temperature substructures on ~ 100 
h^Q kpc scales in the central regions of nine intermediate-redshift (z ~ 0.1) galaxy 
clusters, which resemble those found in the Virgo and Coma Clusters. Each 
substructure contains a clump of hot plasma whose temperature is about 2 — 3 
keV higher than the environment, corresponding to an excess thermal energy of ~ 
]^g58-60 gj,g pgj, ciuxxip. Since if there were no significant non-gravitational heating 
sources, these substructures would have perished in 10®^^ yrs due to thermal 
conduction and turbulent flows, whose velocity is found to range from about 
200 to 400 km s~^, we conclude that the substructures cannot be created and 
sustained by inhomogeneous radiative cooling. We also eliminate the possibilities 
that the temperature substructures are caused by supernova explosions, or by 
the non-thermal X-ray emission due to the inverse-Comptonization of the CMB 
photons. By calculating the rising time of AGN-induced buoyant bubbles, we 
speculate that the intermittent AGN outbursts (> 10^° erg per burst) may have 
played a crucial role in the forming of the high temperature substructures. Our 
results are supported by recent study of McNamara & Nulsen (2007), posing a 
tight observational constraint on future theoretical and numerical studies. 

Subject headings: galaxies: clusters: general — galaxies: intergalactic medium 
— techniques: image processing — X-rays: galaxies: clusters 
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1. INTRODUCTION 

In the frame of hierarchical formation theories, it is beheved that the mass clustering 
develops due to the gravitational instability caused by the intrinsic, tiny Gaussian density 
perturbations that existed in the early universe (see, e.g., Frenk et al. 1996 for a review). 
Since the initial conditions and input physics are assumed to be scale-free, dark halos are 
expected to have evolved self-similarly in time, which can be described by the scaling laws 
over many fundamental parameters (Navarro et al. 1995 and references therein), as have 
been advocated by the numerical simulations in which spherical collapse and adiabatic gas 
dynamics are assumed (e.g., Navarro et al. 1996; Eke et al. 1998). 

On the observational aspect, tight correlations between X-ray luminosity, gas temper- 
ature, gas entropy. X-ray isophotal radius, and virial mass of galaxy clusters have been 
discovered successively with Einstein, ASCA, and ROSAT (e.g., Mushotzky 1984; Mohr & 
Evrard 1997; Horner et al. 1999). These correlations, however, often show notable devia- 
tions from theoretical predictions. For example, in early 1990s, a significant scatter on the 
measured luminosity-temperature curves was reported by, e.g.. Edge & Stewart (1991) and 
Fabian et al. (1994b), who showed that the cool-core galaxy clusters are more luminous at 
a given temperature. The observed mass-temperature relation, on the other hand, shows a 
lower (~ 40%) normalization and a steeper slope with respect to the prediction of adiabatic 
models (e.g., Evrard et al. 1996; Xu et al. 2001). Moreover, an excess in the gas entropy 
was found in the central regions (< 0.05-R2oo) of many low temperature systems (Ponman et 
al. 1999, 2003), which forms a platform of extra energy (~ 1 — 3 keV per particle; Tozzi & 
Norman 2001) on the entropy-temperature curve. These deviations and mismatches, which 
have been confirmed with Chandra and XMM-Newton observations recently (e.g., Arnaud 
et al. 2005; Kotov & Vikhlinin 2006), clearly indicate that our understanding of the thermal 
history of galaxy clusters is far from complete, and the roles played by non-gravitational 
heating processes need to be evaluated properly (e.g.. Kaiser 1991; Evrard & Henry 1991; 
Kauffmann et al. 1993), which will also help solve the coohng flow problem (Fabian 1994a; 
Makishima et al. 2001; Peterson et al. 2003). 

As of today, a number of heating sources have been proposed, which include merger 
shocks, cosmic rays, supernovae, thermal conduction, turbulent dissipation and diffusion, 
and AGN feedbacks (e.g., McNamara & Nulsen 2007; Markevitch & Vikhhnin 2007). How- 
ever, none of the proposed candidates is all-purpose for halting the gas coohng and breaking 

the self-similarities on both galactic and cluster scales. For example, recent observations 
(e.g., Groston et al. 2005; Jetha et al. 2007) showed that, AGN activity, the most represen- 
tative heating source of today, tends to deposit most of its energy into the vast intergalactic 
space, rather than within the host galaxy (< 0.05i?5oo)- Glearly, detailed mappings of gas 
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temperature gradients in the ICM is necessary as a crucial observational constraint in order 
to correctly investigate the gas heating history. 

Two-dimensional gas temperature variations have been observed in many galaxy clusters 
with Chandra and XMM-Newton. By analyzing the temperature map generated with the 
high-quality Chandra data, Fabian et al. (2000) identified a pair of cool gas shells that 
coincide with two radio bubbles, whose radii are ~ 6 and 8 Hjq kpc, respectively, in the 
innermost region of the Perseus Cluster. The authors suggested that the cool gas shells 
were formed as the buoyant gas detached from the terminals of AGN jets and inflated in 
pressure equilibrium with the surrounding gas. Similar phenomenon was also observed in, 
e.g., Centaurus cluster (Fabian et al. 2005), Abell 2052 (Blanton ct al. 2003), and Hydra 
A (Nulsen et al. 2002). Also, in clusters that show clear merger signatures, such as IE 
0657-56 (Govoni et al. 2004), Abell 3667 (Vikhlinin et al. 2001), Abell 665, and Abell 2163 
(Markevitch & Vikhlinin 2001), conspicuous temperature variations were copiously found, 
which are presumably caused by the rapid adiabatic expansion of the shock-heated gas. 

So far, however, in most of the existing works, little has been done to measure the mor- 
phologies of the temperature substructures and to investigate their origins in a quantitative 
way. One of few exceptional works was presented by Shibata et al. (2001), who calculated 
the two point correlation function of the ASCA hardness ratio map of the Virgo Cluster, and 
found that there is a characteristic scale of ~ 300 H^q kpc for the two-dimensional tempera- 
ture variations. The authors associated this characteristic scale with the size of the gas blobs 
arising during the infalling of sub-groups. By analyzing the two-dimensional temperature 
distribution with the XMM-Newton data, Schuecker et al. (2004) reported the presence of 
a turnover at ~ 100 h^^ kpc on the pressure fluctuation spectrum of the Coma Cluster, and 
attributed it to the existence of Kolmogorov/Oboukhov-type turbulent flows. As of today, 
due to the lack of more such quantitative studies, it is not clear if there does exist a preva- 
lence for a characteristic scale ~ 100 H^q kpc in galaxy clusters and groups, and how far it 
can be used to constrain the heating models in the ICM. 

In this work we analyze the Chandra archive data to search for possible characteristic 
scales of temperature variations in a sample of nine intermediate-redshift clusters. In §2 & 3, 

we describe the sample selection criteria and data analysis, respectively. In §4, we discuss the 
physical implications of the results. In §5, we summarize our work. We assume Hq = 70 km 
s~^ Mpc^^, a flat universe for which Qq — 0.3 and Oa — 0.7, and adopt the solar abundance 
standards of Grevesse & Sauval (1998). 
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2. SAMPLE AND DATA PREPARATION 

In order to balance between the angular resolution and detector coverage of the targets, 
we construct our sample with nine intermediate-redshift {z ~ 0.1) galaxy clusters, so that 
the central 400 hjQ kpc regions of the clusters can be fully, or nearly fully covered by the 
S3 CCD of the Chandra advanced CCD imaging spectrometer (ACIS) instrument. All the 
selected clusters have a 0.7 — 8.0 keV flux greater than 5.0 xlO~^^ erg s^^ cm^^. We list 
some basic properties of these clusters in Table 1, which are arranged in the orders of source 
names (col. [1]), names of the cD galaxies (col. [2]), redshifts (col. [3]), richness classes (col. 
[4]), right ascension and declination coordinates (J2000) of the cluster optical centroids (col. 
[5] & [6]), and notes (col. [7]). 

On the optical and infrared images extracted from the (Fig. 1) and 2MAsd!l 

archives, we do not find any strong evidence for pronounced signatures caused by major 
mergers on cluster scales, such as off-center optical and/or infrared substructures. A478, 
A1068, A1650, A2244, and A3112 possess a central galaxy that has been classified as a cD, 
whose size is > 3 times that of any other member galaxy (Rood & Sastry 1971; Takizawa 
et al. 2003). The remaining four clusters are dominated by a pair or a clump of brightest 
member galaxies (Struble & Rood 1987). In addition, Schombert et al. (1989) reported the 
existence of a low- velocity (the radial velocity difference is Av,. ^ 50 km s^^) companion 
galaxy dwelling deep inside the envelope of the cD galaxy of A2244. 

Strong radio emission associated with AGN activity has been detected in four clusters. 
In A478, two symmetric radio lobes are identified at 1.4 GHz, which extend about 3" towards 
northeast and southwest, respectively (Sun et al. 2003). In A1068, two compact radio 
sources are resolved, which are consistent with the locations of the cD galaxy and SDSS 
J104043. 45+395705. 3 (a member galaxy located at about 13" southwest of the cluster center), 
respectively (McNamara et al. 2004). A2204 and A3 112 are reported to possess an extended 
radio halo that roughly traces the spatial distribution of the X-ray gas, respectively (Sanders 
et al. 2005; Takizawa et al. 2003). In the remaining five clusters, neither radio lobes nor 
luminous radio halos has been found in the central region in previous works, suggesting that 
they have not experienced strong AGN outbursts during past few tens of Myr, if the life 
time of such radio sources can be estimated by tgync = ^^b^^ ^^"^ -B ~ 10 /xG (Takahashi 
& Yamashita 2003; Donahue et al. 2005). As a unique case in our sample, A2556 harbors 
an off-center extended radio source at about 40" northeast of the center, which is often 



^http:/ /archive. stsci.edu/dss/ 

^ ht tp : / / www . ipac. caltech . edu / 2mass / 



- 5 - 



speculated to be the relic of a merger-driven shock (e.g., Govoni et al. 2004). 

All the X-ray data analyzed in this work were acquired with the ACIS S3 CCD (Table 
2), with the focal plane temperature set to be -120 °C. Using the CIAO v4.0 software 
and CALDB v3.5.0, we removed the bad pixels and columns, as well as events with ASCA 
grades 1,5, and 7, and then executed the gain, CTI, and astrometry corrections. In order to 
identify possible strong background flares, lightcurves were extracted from regions sufficiently 
far away, i.e., > 4' from the X-ray peaks. Time intervals during which the count rate exceeds 
the average quiescent value by 20 percent were excluded. The SI chip data were used to 
crosscheck the results when it was available. 



3. X-RAY IMAGING AND SPECTRAL ANALYSIS 

3.1. X-Ray Images and Surface Brightness Profiles 

Figure 1 shows the intensity contours of the X-ray emission in 0.3 — 10 keV of our 
sample clusters, which overlay the corresponding optical DSS images. In all cases, the X- 
ray morphology appears to be smooth and symmetric on > 100 h^Q kpc scales, showing no 
remarkable twists, fronts, or edges that indicate recent major mergers. On smaller scales, 
however, there exist substructures, such as X-ray cavities coinciding with the AGN lobes 
(Sun et al. 2003), stripes, filaments, and arc-like sharp fronts (Wise et al. 2004; Sanders et 
al. 2005). 

For each cluster, we examine the vignetting-corrected surface brightness profile (SBP) 
extracted in 0.7 — 8.0 keV, after filtering all the point sources that are detected above the 
3(T level of the local background by using both the cclldetcct and wavdctect tools with their 
default settings (Fig. 2). By applying the x^-test, we find that, except for the SBP of 
A2204, which can be well described with a single-/? model S{r) = S'(0)(1 + (r/r^)^)^'^'^^^/^, 
all the obtained SBPs are consistent with the empirical two-/3 model S{r) = 5'i(0)(l -|- 
(r/rci)2)-^'^i+^/2 + S2{Q){1 + (r/rc2)^)"^^^+^/^ (e.g., Ikebe et al. 1996), because in the central 
10-20 /ifo^ kpc (A1650, A2556, and A3112) or 20-40 h^^ kpc (A478, A1068, A1201, A2104, 
and A2244) there exists an emission excess beyond the /3 model that can best describe the 
SBPs of outer regions. Such an excess is often seen in clusters that host a central dominating 
galaxy, which is usually classified as a cD (Makishima et al. 2001). 
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3.2. Azimuthally-averaged Spectral Analysis 

3.2.1. Background 

In the spectral analysis that follows, we utilized the Chandra blank-sky templates for 
the S3 CCD as the background. The templates were tailored to match the actual pointings, 
and the background spectra were extracted and processed identically to the source spectra 
(§3.2.2). After rescaling each background spectrum by normalizing its high energy end to 
the corresponding observed spectrum, we further enhance the accuracy of the background by 
adding an unabsorbed APEC model with fixed temperature (0.2 keV) and metal abundance 
(1 Zq] Snowden et al. 1998) to account for variations in the Galactic X-ray background 
(GXB; Vikhlinin et al. 2005). We note that the GXB component is significant only in 
A2204, where it has a flux of 2.0 xlO~^^ erg s"^ cm~^ arcmin"^ in 0.3 — 1.0 keV and is 
nearly uniform over the S3 chip. This was also reported by Sanders et al. (2009), who 
indicated that it is probably due to the excess Galactic foreground. 

3.2.2. Models and Results 

Here we calculate the azimuthally averaged temperature profiles, which will be used to 
create the reference temperature maps applied in wavelet detection (§3.3.2) and calculation 
of excess energies in the temperature substructures (§4). We divided the inner 400 h^^ kpc of 
each cluster into concentric annuli, masked all the detected point sources (~ 15 point sources 
per cluster), and extracted the spectrum from each annulus. The corresponding spectral 
redistribution matrix files (RMFs) and auxiliary response files (ARFs) were generated using 
the ClAO tool mkwarf and mkacisrmf. The spectra were grouped with > 20 counts per 
bin to allow the use of statistics. We used XSPEC 12.4.0 package to fit the spectra, by 
modeling the gas emission with an absorbed APEC component and employing the PROJCT 
model to correct the projection effect. The absorption column density Ah was fixed to the 
Galactic value (Dickey & Lockman 1990) in all cases, except for A478, for which Ah was set 
to be about twice the Galactic value (Sanderson et al. 2005; Table 3); since the absorption 
excess in A478 exhibits a spatial extension beyond the cluster's core region, Pointecouteau 
et al. (2004) argued that it is likely to be associated with a local absorption substructure in 
our Galaxy. In order to calculate the error ranges of the model parameters, we executed the 
XSPEC command steppar and carried out multiple iterations at every step, to ensure that 
the actual minimum is found. The best-fit gas temperature and abundance profiles are 
illustrated in Figure 3, in which multiple sets of annuli are used for each cluster to crosscheck 
the results. We note that the inclusion of an additional gas component could only slightly 
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improve the fits for the central regions of A478, A1650, A2204, A2244, and A3112, which is 
insignificant as indicated by the F-test. 

A significant temperature dechne by a factor of > 2 is observed in the central regions 
of A478, A1068, A2204, A2556, and A3112, four of which (except A2556) possess radio 
halos and/or lobes in the center (§2). In the remaining four clusters without apparent radio 
emission, the inward temperature drop is somewhat ambiguous, as the difference between the 
core temperature and mean cluster temperature, which refer to the temperatures measured 
in the central 100 hjQ kpc (~ O.lrsoo) and 100 — 400 hjQ kpc regions (Fig. 3), respectively, 
is found within Scr error ranges. If we apply the classification criteria of Sanderson et al. 
(2006), the latter four clusters are non-cooling-core clusters, whose central regions have not 
yet fully condensed. These results refiect the trend for the cooling core clusters to host 
central radio sources (Burns 1990). As shown in Figure 3, we also adopt Eqs.(4) — (6) of 
Vikhlinin et al. (2006) to describe the obtained radial temperature profiles in a smooth 
form, which will be used in §3.3.2 to calculate the reference temperature maps in order to 
examine the significances of detected temperature substructures, and in §4 to calculate the 
excess energies in these substructures. 

In A478, A1068, A1650, A2204, A2556, and A3112, the abundance increases inwards 
significantly (68% confidence level) from ^ 0.2 Zq in the outermost annuli to > 0.6 Zq in 
the innermost regions (< 30 h^Q kpc). Emission measure- weighted gas temperature and 
abundance of each cluster are also calculated (Table 3), which are in good agreement with 
previous results obtained with Chandra and XMM-Newton (Sun et al. 2003; Takahashi & 
Yamashita 2003; Takizawa et al. 2003; Pointecouteau et al. 2004; Wise et al. 2004; Donahue 
et al. 2005; Sanders et al. 2005; Sanderson et al. 2005). 

3.3. Two-Dimensional Analysis 

3.3.1. Projected Temperature Maps 

Following the methods of, e.g., O'SuUivan et al. (2005), Maughan et al. (2006), and 
Kanov et al. (2006), we define a set of knots (> 5000 knots per cluster) in the central 
300 /ifo^ kpc region, which are randomly distributed with a separation of Ajj < 5" between 
any two adjacent knots i and j. To each knot we assign an adaptively sized circular cell, 
which is centered on the knot and contains more than 800 counts in 0.7 — 8 keV after the 
background is subtracted. The typical cell radius ranges from < 5" at the center to about 
10" at r > 150 h^^ kpc, which is always larger than the associated Aj^. We calculate 
the projected temperature of each cell Tc(ri), which is assigned to the corresponding knot. 
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by fitting the spectrum extracted from the cell with an absorbed APEC model, using the 
absorption column given in Table 3; adding another thermal component can neither improve 
the fit significantly, nor make the parameters constrained more tightly. In the fittings, the 
source spectrum is grouped to have > 20 photons in each energy bin. In addition to the 
standard statistics, we also have applied the maximum likelihood statistics, and find that 
the best-fit parameters obtained with the two statistics agree nicely with each other within 
the la error ranges. The calculated la errors of gas temperature are typically < 10% in 
the central 100 h^^ kpc, which increase progressively to < 15% (A478, A1650, A2244, and 
A3112) or < 20% (the remaining five clusters) in 100 - 300 /i^ kpc. 

For any given position r, we define a scale s(r), so that there are at least 800 counts 
contained in a circular region centered at r, whose radius is s(r). We calculate the projected 
temperature at r by integrating over all the knots ri located in the circular region, 

r(r) = 5^(Gr,(i?r,rjTe(rO)/5^Gr,(i?r,rJ, whcu i^^.r. < ^(r), (1) 
rj ri 

where Rr^n is the distance from r to ri, and Gr^ is the Gaussian kernel whose scale parameter 
a is fixed at s(ri). Since s(r) is essentially proportional to the square root of the local counts, 
the obtained temperature map T(r) (Fig. 4) is less affected by the statistical uncertainties 
caused by surface brightness fluctuations. In addition, angular resolutions of ~ 10" are 
guaranteed by the use of compact Gaussian kernel. Maps for the upper and lower limits of 
T{v) are calculated in a similar way. 



3.3.2. Temperature Substructures 

Wavelet Transform 

A visual inspection of the obtained temperature maps (Fig. 4) indicates that there exist 
significant substructures, which typically have linear sizes of ~ 100 K^q kpc, in all of the 
sample clusters. In order to describe these substructures in a quantitative way, we analyze 
the statistical properties of the temperature maps by applying the wavelet transform, a 
powerful tool in imaging process that identifies and disentangles small embedded features 
as a function of scales (e.g., Vikhlinin et al. 1997; Starck & Pierre 1998). In an analogue 
study, Bourdin & Mazzotta (2008) conducted the wavelet transform to achieve an optimized 
spatial binning of temperature maps. We employ a non-orthogonal version of discrete wavelet 
transform, the a Trous wavelet (Shensa 1992; Starck et al. 1995; Vikhlinin et al. 1998), which 
defines a straightforward decomposition rule. To be specific, at position r = (,t, y) and scale 
level i {i — 1,2,...), the a Trous wavelet Wi{x,y) equals the signal difference between two 
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smoothed maps Ci_i(x, y) — Ci{x, y), where co{x, y) is the obtained temperature map T(x, y), 
and Ci{x,y) connects with Ci_i(a;, |/) via Ci{x,y) = ^^^h{m,n)ci-i{x + 2^~^m,y + 2^~^n), in 
which /i(m, n) is the convolution kernel, and (m, n) is its grid. Here /i(m, n) is chosen to be 
a compact Gaussian h{m,n) = e~(™^+"'^)/^/2'^. Each temperature map is thus decomposed 
into a set of subimages Wi{x, y) on nine scales (2* pixels, i = 1, 2, ..9), respectively. 

Based on the obtained Wi{x,y) maps, we calculate the significant wavelet coefficient 
Wi{x,y) to detect the temperature substructures on different scales. For each cluster, we 
create a set of random knots ri as used in §3.3.1, and assign a randomized temperature T^{ri) 
to each knot, which is constrained by the measured temperature profile T(r) and its error 
ranges (Fig. 3). By applying Eq.(l) to the T^'('^i) iiiaP) we obtain a reference temperature 
map Ti.ef(r), in which the effects of the radial temperature gradients and knot tessellation 
are both involved. A series of such reference temperature maps are created by repeating this 
random process, and each reference map is decomposed into subimages Wi^refi^, y) on nine 
scales by performing the same a Trous wavelet transform as described above. On each scale 
i, we determine the significant coefficient as Wi{x,y) = Wi{x,y) when Wi{x,y) > 3o"j,rcfj and 
Wi{x,y) = when Wi{x,y) < Sai^^ei (e.g., Slezak et al. 1994; Starck et al. 2003), where 
(Tj i-cf is the standard deviation of the simulated Wi^^ci{x,y) maps so that 3(Tj^rcf is equivalent 
to a confidence interval of 1.3 x 10^^ for Gaussian noise. A similar method was used to 
determine the actual significances of diffuse sources on ROSAT images (e.g., Rosati et al. 
1995; Grebenev et al. 1995; Damiani et al. 1997). 

Since the a Trous transform is essentially a non-orthogonal algorithm, the obtained 
coefficient Wi{x,y), and thus Wi{x,y), may be systematically biased from the exact solution 
(Starck & Pierre 1998). This bias can be corrected by using an iterative method (Van Cittert 
1931; Biviano et al. 1996; Starck et al. 2003; Bourdin et al. 2004), the iteration routine of 
which is defined as 

where W"(a;, y) is the corrected detection on scale i after n iteration(s) {W^{x^ y) = Wi{x, y) 
is the initial run), and P is a non- linear synthesized operator of inverse wavelet transform, 
wavelet transform, and filtering operation using 3(7j,ref as the threshold. Typically the calcu- 
lation converges within five iterations, which yields a corrected Wj(x, y) map as the unbiased 
reconstruction of the temperature substructures on scale i (Fig. 4). We define the wavelet 
spectrum as (| Wi{x,y) p) (Torrence & Compo 1998), and plot it in Figure 5 as a function 
of scale for each cluster. 



As a crosscheck, we also employ the two-dimensional B-spline function (e.g., Bourdin 
& Mazzotta 2008) as the scaling function 0(x, y) to perform a compact a Trous trans- 
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form in the Fourier domain {u,v), so that the Fourier transform of the wavelet coefficient 
Wi{x,y) is directly calculated as Wi{u,v) = (1 — h{2^~^u,2^~^v))ci-i{u,v), where the filter 
h{u,v) — (f){2u,2v)/4){u,v) when \u\, \v\ < 1/2 (in this paper the Fourier transform of any 
function / is written as /). By carrying out the inverse Fourier transform to Wi{u, v), Wi{x, y) 
maps are obtained, and then used to calculate the corrected Wi{x, y) maps following the pro- 
cedure described above. The corresponding wavelet spectra are illustrated in Figure 5 for 
comparison, which agree nicely with those obtained with the Gaussian wavelet. We find that 
the most popular scales of the projected temperature substructures are in the range of about 
50 — 200 h^Q kpc, and about 70% of these substructures are located at 100 - 200 K^q kpc 
from the cluster center, exhibiting irregular shapes and spatial distributions (Fig. 4). 

Power Spectrum 

As a straightforward approach, we have also studied the power spectra of the two- 
dimensional gas temperature distributions. Adopting the method of Schuecker et al. (2004), 
we determine the flat flelds T'(r) by applying a low-pass wavelet fllter with a scale of 2^ pixels 
on the original temperature maps, and calculate the normalized temperature fluctuation 
distributions as ^^(r) = T{r)/T{r) — 1. By carrying out the Fourier transform of ST{r) 

5(k) = l^(5r(r)e-''dr, (3) 

where S is the projected area, we obtain the power spectrum P{k), 

P{k) = {\ 6{k) \') . (4) 

Following the method of Schuecker ct al. (2004), we carry out Monte-Carlo simulations 
for each cluster to estimate the temperature fluctuation at any given knot ri based on the 
temperature error range measured at the same knot (§3.3.1). After a randomized value 
Terr(ri) is assigucd to each knot, a random temperature error map Terr(r) is created by 
applying Eq.(l) to the Terr(ri) map, in which the uncertainties introduced by data statistics, 
spectral model fltting, and cell tessellation are all involved. Then, by combining T(r) and 
Te]-r(r) maps we create a random temperature map Tj.g^j^{r). The la error bars shown in 
Figure 6 are calculated from the variances of power spectra obtained from ten such Tran(r) 
maps for each cluster. In general, the power spectra P{k) show a signiflcant peak in 50 — 200 
hjQ kpc, and drop quickly beyond 200 — 250 h^Q kpc. This strongly supports our conclusion 
derived with the wavelet analysis (Fig. 5). 
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4. DISCUSSION 

We present robust evidence for the prevailing existence of hot gas clumps, whose char- 
acteristic scales are ~ 100 h^^ kpc, in the central regions of a sample of nine intermediate- 
redshift galaxy clusters. These substructures are not caused by uncertainties in background 
subtraction, response calibration, and spectral model fittings, nor are they related to back- 
ground flares, the effects of which have been largely corrected and/or flxed. The possibihty 
of the temperature substructures arising from the non-thermal emission of low-mass X-ray 
binaries (LMXBs) can also be eliminated, since the ICM is only sparsely filled with galaxies. 
Our results confirm the previous findings in the Coma Cluster (Schuecker et al. 2004) and 
the Virgo Cluster (Shibata et al. 2001). 

We estimate the excess thermal energy E'excess in each high temperature substructure 
(Table 4), where the gas temperature is typically higher than environment by ATx ~ 2 — 3 
keV (Fig. 4), by defining the thermal energy as the product of gas pressure P and volume 
V, and assuming a spherical or elliptical geometry for the volume of the hot gas clump. We 
obtain -^excess = X V = ATxTicV ~ X0^^~^° erg for each clump, where ric is the electron 
density determined in the deprojected spectral analysis (§3.2.2). To account for azimuthal 
variations in the electron density, we include a 10% systematic error (Sanders et al. 2004) 
in the calculations. As the projection effect is not corrected for ATx, the obtained excess 
energy may have been underestimated by a factor of about two. The origin of this energy 
excess is discussed as follows. 



4.1. Inhomogeneous Cooling, Conduction, and Turbulence 

The detected temperature substructures may have been formed due to the spatially 
inhomogeneous cooling in the ICM. Since for a unit volume, the radiative cooling rate is 
-Rrad — TT'e^-^radiT, Z), where Arad IS the cooliug fuuction calculated by including both con- 
tinuum and line emissions, the radiative cooling time of the substructures can be estimated 
as 

Typically, the calculated coohng time is ~ 10^~^° yrs, and is thus nearly an order of mag- 
nitude longer than the characteristic time for thermal conduction alone to smear out the 
temperature substructures (Table 4), which is given by 



5h . ^neksSTx 
— 5r = ^ 



^cond ^ —Sr = 4 d7V (6) 
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where 5h is the enthalpy excess per unit volume in the hot substructure, q is the conductive 
flux, 5r is the characteristic scale length of the substructure, / = 0.2 is the factor to assess 
the suppression caused by magnetic fields (e.g., Narayan & Medvedev 2001), and /«cond is the 
thermal conductivity of hydrogen plasma. 



Keond ^ 5.0 X 10-^t| [erg cm-^ s'^ K'^] (7) 

for Coulomb gas (Spitzer 1962). This indicates that, unless the conduction process is ex- 
tremely hindered (i.e., when / ^ 1; Maron et al. 2004; Markevitch et al. 2003), the effect 
of the radiative cooling is less important in the formation and evolution of temperature 
substructures. Hence the possibility of the substructures being formed by inhomogeneous 
cooling can be excluded. 

The destruction of a temperature substructure via conduction can be accelerated by tur- 
bulence (Dennis & Chandran 2005), whose effect can be evaluated by introducing a modified 
conductivity (Cho et al. 2003; Voigt & Fabian 2004) 



1/2 

<ond = anMl ( ) , (8) 



3//mp 

where a is the ratio of the turbulence velocity Vturb to the local adiabatic sound speed Cg 
{a — Vturb/cs)) a-nd I is the turbulence length scale. To calculate k'^^j^^, the turbulence velocity 
Vtmh can be estimated as follows. As there is evidence that heating continuously compensates 
a large fraction of radiative coohng since z~ 0.4 (Bauer et al. 2005; McNamara & Nulsen 
2007), we may assume an approximate energy balance between heating and cooling for the 
cluster (e.g., Kim & Narayan 2003; Zakamska & Narayan 2003; Dennis & Chandran 2005) 
as 

Rc^H + T + Q, (9) 

where R is the radiative loss, and H, T and Q are the heating rates of thermal conduction, 
viscous dissipation and turbulent diffusion, respectively. Since 

= V-(«cond/Vr), (10) 

3 

p _ CdissP'^turb ^-|^-|^'^ 

and 

V- (DpTVs), (12) 

where p is the gas mass density, s is the specific entropy, Cdiss ~ 0.42 is a dimensionless 
constant (Dennis & Chandran 2005 and references therein), and D 10^^ cm s~^ is the 
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diffusion coefficient (Rebusco et al. 2006). By subtracting Eqs.(lO) — (12) into Eq.(9), we 
fiave 

dr r dr r CdissP 

if we define \1/ = //tcond^ and = (Dennis & Cliandran 2005). Using tlie observed 

deprojected gas density and temperature profiles (Fig. 3), we calculate Vtmh and show the 
results as a function of turbulence scale I in Figure 7. If we follow Schuecker et al. (2004) 
to adopt the lengths indicated by the turnover points on the wavelet spectra (Fig. 5) as the 
turbulence scale, we find that Vtmh commonly ranges from 200 to 400 km s~^, or a ~ 0.2 — 0.4, 
for our sample. By using Eqs.(6) and (8), we find that the inclusion of turbulence reduces 
the characteristic time tcond to the modified characteristic time t'^^^^^ by a factor of about 3 
(Table 4), which sets an upper limit (~ 10® yrs) on the duty cycles of the dominating heating 
sources. 



4.2. AGN Feedback 

The calculated i?cxcess (Table 4) suggests that the hot clumps may be dying or remnants 
of buoyant bubbles, into which the central AGN has injected energy via shocks and turbu- 
lences (e.g., McNamara & Nulsen 2007). Such AGN- induced hot gas clumps are predicted 
in numerical simulations (e.g., Quilis et al. 2001; Dalla Vecchia et al. 2004; Vernaleo & 
Reynolds 2006) and theoretical calculations (e.g., Pizzolato & Soker 2006), whose scales are 
expected to be about 100 kpc by the end of the buoyancy process. By examining the tem- 
perature maps (Fig. 4), we note that some hot clumps appear in pairs on both sides of the 
central galaxies (#1 and #2 in A1201; #1 and #2 in A2204; #1 and #3 in A3112), which 
suggests possible AGN activity taking place at the cluster cores. However, can an AGN- 
heated temperature substructure survive conduction and turbulence destructions? To answer 
this question, we estimate the characteristic rising time of a bubble as tbuoy = ''"b/'^buoyi where 
Tb is the travel distance measured from the cluster center to the location of a hot clump, 
and Vbuoy is the travel velocity. In terms of the Keplerian velocity t'K('"b), a large bubble is 
expected to travel through the ambient gas at 



^^buoy(rb) ^ \/^^^K(rb), (14) 

where C ~ 0.75 is the drag coefficient and it!b is the bubble radius in kpc (Ghurazov et al. 

2001). This yields tbuoy smaller or very close to t^ond ^'^^ ~ 95% of the substructures (Table 
4), indicating that such substructures can survive in most cases. Even for the clumps with 
^buoy > ^cond' ^hc AGN heating scenario may still stand, if the distance traveled by the bubble 
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is shorter (e.g., Gardini 2007). It should be noted that no hot gas clump is found associated 
with any known radio source in all sample clusters. This, however, does not contradict with 
the AGN heating scenario, because the lack of radio sources can be explained by the fact 
that the relativistic electrons may have lost most of their energy when the bubbles arose, as 
inferred by t^ync < ^buoy (Wnc ~ 10'^"^ yrs; §2). 

4.3. Other Possible Mechanisms 

Numerical and theoretical works show that shocks caused by supersonic galaxy motion 
in the ICM can generate observable temperature substructures via gas compression and 
friction (e.g., El-Zant ct al. 2004), and even the low-speed galaxy infalling may also be a 
heating source candidate due to the efficient energy conversion via the magnet ohydro dynamic 
turbulence and magnetic rcconncction (Makishima et al. 2001; reviewed by Markevitch & 
Vikhlinin 2007). The energy dissipation rate in the central region of a cluster is as large as 
10^^ erg s~^ (Fujita et al. 2004), which assembles ~ 10^^ erg during the merger process and 
is thus sufficient to account for the excess energy in the observed temperature substructures. 
Nevertheless, no convincing evidence is given in numerical studies that the non-filamentary 
morphologies of the observed temperature substructures can appear commonly in merger 
process (e.g., Poole et al. 2006). 

We also have estimated the contribution of supernova feedback to the IGM heating 
over the substructure lifetime {t'^^^^^) using the observed rate of type la supernovae (SNe la; 
Sharon et al. 2007), which implies that about 0.8 — 6.1 x 10"^ SN la explosions have occurred 
in the central 400 hjQ kpc region of each sample cluster. Since the averaged dynamic energy 
injection into the environment is 4 x 10''"''^ erg per SN la event (Spitzer 1978), and ~ 10 % of 
this energy is assumed to have been used to heat the gas (e.g., Thornton et al. 1998), the 
total SN la heating is found to be < 3 x 10^^ erg per cluster per 10^ yrs. The contribution 
of type II supernovae (SNe II) will not significantly increase the total supernova heating, 
since even in the IR-luminous cluster A1068, where the ratio of SN II rate to SN la rate is 
3 (Wise et al. 2004), the SN II heating is still lower than 8 x 10^'^ erg. Therefore, we 
conclude that the supernova heating is insufficient to create the observed high temperature 
substructures. 

High temperature substructures may also occur as a result of the inverse-Gomptonization 
(IG) of the cosmic microwave background (CMB) photons that are scattered by the rela- 
tivistic electrons in the IGM, since in a Chandra ACIS spectrum it is difficult to distinguish 
between such an IG component and a high temperature thermal component (Pctrosian et 
al. 2008). If in the IGM 1% of the electrons are relativistic (7 - 10^; Eilek 2003), we will 



- 15 - 



have an IC energy conversion rate of 6(7) = |^7^C^cmb — 6 x 10^^ erg s"^, where iTt is the 
Thomson cross section and Uqub is the CMB energy density at the location of the cluster 
(Sarazin 1999), which implies that the flux of the IC component is sufficiently high to bias 
the temperature measurements by about 1 — 3 keV. However, a tight correlation between the 
high temperature substructures and radio synchrotron sources is required in this scenario, 
which is actually not observed, possibly due to the lack of relativistic electrons. Therefore, 
the IC mechanism for the formation of the temperature substructures can be eliminated. 



5. SUMMARY 

We reveal the prevailing existence of temperature substructures on ~ 100 hjt^ kpc 
scales in the central regions of nine intcrmcdiatc-rcdshift galaxy clusters. By comparing the 
characteristic destruction times of the temperature substructures with the gas cooling times 
and the rising times of AGN-induced bubbles, we conclude that the AGN outbursts may 
have played a crucial role in the forming of the observed temperature substructures. Our 
results agree with those found earher in the Virgo and Coma Clusters (Shibata et al. 2001; 
Schuecker et al. 2004). 
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Table 1. Galaxy Clusters in Our Sample 



Name 


cD Galaxy 


Redshift 


R» 


RA'' 
(h d s; J2000) 


DEC'' 
(d m s; J2000) 


Notes^ 


A0478 


2MASX J04132526+1027551 


0.0881 


2 


04 13 20.7 


+ 10 28 35 


1,3 


A1068 


2MASX J10404446+3957117 


0.1375 


1 


10 40 47.1 


+39 57 19 


2,3 


A1201 




0.1688 


2 


11 13 01.1 


+13 25 40 




A1650 


2MASX J12584149-0145410 


0.0845 


2 


12 58 46.2 


-01 45 11 




A2104 




0.1554 


2 


15 40 06.8 


-03 17 39 




A2204 




0.1523 


3 


16 32 45.7 


+05 34 43 


2,3 


A2244 


FIRST J170242.5+340337 


0.0968 


2 


17 02 44.0 


+34 02 48 




A2556 




0.0871 


1 


23 13 01.6 


-21 37 59 


4 


A3112 


ESQ 248- G 006 


0.0750 


2 


03 17 52.4 


-44 14 35 


3 



"Abell richness class (Abell et al. 1989). 
''Positions of the cluster optical centroids. 

'^X-ray and radio substructures in the clusters, including 1 — X-ray cavities, 2 — cold fronts, 3 — central 
radio sources, and 4 — off-center radio source. 
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Table 2. Observation Log 



Target 


ObsID 


ACIS CCD 


Mode 


Date 


Raw/Clean Exposure 










dd mm yyyy 


(ks) 


A0478 


1669 


235678 


FAINT 


27/01/2001 


42.9/41.9 


A1068 


1652 


236789 


FAINT 


04/02/2001 


27.2/26.0 


A1201 


4216 


35678 


VFAINT 


01/11/2003 


40.2/25.4 


A1650 


4178 


35678 


VFAINT 


03/08/2003 


27.6/26.0 


A2104 


895 


235678 


FAINT 


25/05/2000 


49.8/48.1 


A2204 


499 


235678 


FAINT 


29/07/2000 


10.2/9.3 


A2244 


4179 


35678 


VFAINT 


10/10/2003 


57.7/56.0 


A2556 


2226 


35678 


VFAINT 


05/10/2001 


20.2/19.8 


A3112 


2516 


35678 


VFAINT 


15/09/2001 


17.2/15.1 
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Table 3. Global Properties of the Halo Gas' 



Name 


Aperture 


kT 


Abundance 








(keV) 


(^0) 


(1020 cm-2) 


A0478 


405.4 




0.47i°;°i 


26.75to 5^ 


A1068 


418.7 


'-'•"'-'-0.10 


56+°-°'' 


1.40 


A1201 


425.5 


4-8ll°o:g 


S2+°-" 
'-'••^^-0.09 


1.61 


A1650 


390.5 


5.96l°Jt 


'J-48_o.08 


1.56 


A2104 


397.6 


y.bu_o 7g 


n 47+0-07 

-0.08 


8.69 


A2204 


391.0 


7 24+0-*3 


51+0-06 


5.67 


A2244 


396.9 


C r,c-|-0.10 

J.^J-0.08 


„ 20+0-03 
'J-^y-0.03 


2.14 


A2556 


401.3 


'^•-'■'-0.10 


n 4S+0-05 

^•^■5-0. 04 


2.05 


A3112 


385.5 


. QO-H0.13 
4t.UJ_Q ^2 


n. 7Q-I-0.08 

'-'.19-0.07 


2.61 



*An absorbed APEC model is used to estimate the emission 
measure-weighted gas temperatures kT and abundances. The error 
bars axe given at 90% confidence level. 

''Column densities are fixed to the Galaxitic values (Dickey & Lock- 
man 1990), except for A478, which shows an absorption larger than 
the Galactic value (15.1 xlO^o cm-2). 



-24- 



Table 4. Characteristic Timescales for the Detected Temperature Substructures 



No." 


r ^ 




-I— 'excsss 


tcool 


^cond 


t' s 

cond 


t h 

^buoy 










('10® vrs^ 


do® vrs^ 


flO* \TS) 


do® vrs^ 


A0478[l] 


95 


0.016 


1.9 ±0.8 


16.4 ± 5.2 


7.7 ±2.0 


1.9 ±0.4 


1.6 


[2] 


141 


0.010 


7.4 ± 2.2 


23.8 ±6.1 


3.2 ± 0.8 


1.7±0.6 


2.0 


[3] 


155 


0.009 


7.2 ±3.7 


30.0 ±7.4 


2.4 ±0.4 


1.7 ±0.6 


2.2 


[4] 


157 


0.009 


4.3 ±2.5 


33.6 ±8.1 


4.7 ±0.8 


2.3 ±0.5 


2.3 


A1068[l] 


184 


0.004 


11.0 ±4.8 


54.6 ± 13.8 


16.7 ±2.2 


2.2 ±0.8 


2.7 


[2] 


237 


0.002 


14.6 ±8.7 


106.2 ± 19.8 


8.7 ±2.2 


5.4 ±1.5 


3.2 


[3] 


244 


0.002 


32.2 ± 15.8 


135.9 ± 30.4 


10.3 ±2.9 


5.3 ± 1.6 


3.3 


A1201[l] 


123 


0.005 


17.4 ± 10.8 


69.1 ± 16.7 


5.8 ± 1.2 


4.2 ±0.9 


2.0 


[2] 


171 


0.003 


36.0 ± 20.7 


87.6 ± 17.8 


11.1 ±2.1 


5.9 ±1.4 


2.7 


[3] 


185 


0.003 


29.3 ± 16.7 


118.4 ± 30.5 


9.2 ±1.7 


4.7 ±1.0 


2.9 


A1650[l] 


72 


0.013 


4.6 ±2.0 


26.1 ± 12.0 


3.2 ±1.2 


1.5 ±0.8 


2.4 


[2] 


108 


0.007 


12.2 ±5.2 


42.3 ± 12.9 


6.3 ±2.3 


2.2 ±0.5 


2.9 


[3] 


113 


0.007 


4.6 ± 1.5 


38.7 ± 12.2 


4.7 ± 1.4 


2.3 ±0.5 


3.0 


[4] 


171 


0.005 


12.6 ±4.5 


61.5 ± 13.6 


2.8 ± 1.1 


2.5 ±1.1 


3.6 


A2104[l] 


132 


0.007 


7.0 ±5.6 


53.7 ± 16.8 


3.4 ±0.7 


2.7 ±0.5 


2.8 


[2] 


253 


0.003 


47.8 ± 26.5 


87.7 ± 17.8 


1.8 ±0.5 


0.4 ± 0.3 


4.0 


[3] 


261 


0.003 


41.5 ±23.9 


122.8 ±31.2 


1.7 ±0.5 


0.2 ±0.3 


4.1 


A2204[l] 


94 


0.017 


65.8 ±26.1 


22.8 ±6.5 


1.5 ± 0.5 


1.4 ±0.5 


1.4 


[2] 


117 


0.013 


54.0 ±26.8 


22.8 ±6.8 


3.5 ± 0.6 


1.9 ±0.5 


1.6 


[3] 


172 


0.009 


145.7 ± 72.3 


47.2 ± 11.7 


6.0 ± 1.6 


2.9 ±1.0 


2.0 


[4] 


220 


0.006 


104.4 ±42.1 


61.8 ± 13.9 


2.8 ±0.9 


2.1 ±0.8 


2.4 


A2244[l] 


74 


0.014 


3.6 ± 1.0 


22.7 ±6.1 


15.0 ±0.9 


2.5 ±0.5 


1.3 


[2] 


100 


0.010 


8.5 ±3.2 


29.8 ± 7.8 


8.9 ±0.7 


2.5 ±0.5 


1.7 


[3] 


155 


0.006 


4.9 ±2.1 


42.0 ± 13.0 


6.4 ±0.8 


2.5 ±0.6 


2.3 


A2556[l] 


108 


0.005 


13.6 ±3.8 


62.3 ± 12.8 


13.9 ±2.1 


3.7 ±0.8 


1.8 


[2] 


141 


0.004 


13.0 ±2.2 


70.8 ± 12.0 


10.7 ± 1.5 


3.2 ±0.7 


2.3 


[3] 


144 


0.004 


7.0 ± 3.6 


65.3 ± 15.4 


16.0 ±2.6 


3.6 ±0.9 


2.3 


A3112[l] 


73 


0.014 


2.9 ± 1.3 


20.2 ± 6.4 


9.5 ±1.1 


1.9 ±0.6 


1.2 


[2] 


121 


0.006 


5.2 ±2.2 


45.1 ±9.1 


11.4 ± 1.5 


2.9 ±0.7 


1.8 


[3] 


125 


0.006 


6.8 ±3.6 


41.9 ±9.5 


7.2 ± 1.2 


2.2 ±0.6 


1.9 


[4] 


141 


0.006 


3.8 ±2.2 


39.6 ±9.6 


2.8 ±0.8 


2.1 ±0.6 


2.0 



'^Detected temperature substructures. See Figure 4. 

''Distances between the geometric centers of the substructures and the cluster centers. 
"^Averaged electron densities of the substructures. 
"^Excess thermal energies contained in the substructures. 
''Radiative cooling times of the substructures. 

'Destruction times of the substructures due to thermal conduction (Eq.(6)). 

^Modified destruction times of the substructures due to thermal conduction, viscous dissipation, and turbu- 
lent diffusion (Eq.(8)). 

''Rising times of the temperature substructures if they are AGN-induced bubbles. 
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Fig. 1. — DSS optical images of the sample clusters, on which the 0.3 — 10.0 keV X-ray 
intensity contours are plotted. To calculate the X-ray contours we subtract the background, 
correct the exposure, and smooth the image with a minimum signal-to- noise (S/N) of 3 and 
a maximum S/N of 5 per smoothing beam. 
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Fig. 2. — Radial surface brightness profiles extracted in 0.7 — 8.0 keV, which are corrected 
for both background and exposure. The best-fit two-/? and f3 models are shown with solid 
and dash lines, respectively. For A2204 only the j3 model is plotted, which is sufficient to 
describe its SBP. 
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Fig. 3. — Projected and deprojected temperature {upper panel) and abundance profiles 
{lower panel) for two sets of annuli. Error bars are given at 90% and 68% confidence levels 
for temperature and abundance profiles, respectively. 
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R, A. (J2000) 

Fig. 4. — Projected temperature maps, on which the most significant substructures with the 
characteristic scales indicated by the wavelet spectra (Fig. 5), i.e., i = 7 (100 — 200 /ifo^ kpc; 
solid line) and i = 6 (50 — 100 hj^ kpc; dotted line) are marked (§3.3.2). The representative 
substructures are labeled with the same numbers as used in Table 4. 
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Fig. 5. — Wavelet spectra of the two-dimensional temperature distributions T{r) derived 
with both the Gaussian wavelet (cross) and the B-spline wavelet (diamond). Error bars are 
given at la confidence level (§3.3.2). 
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Fig. 6. — Power spectra of the two-dimensional temperature distributions T(r) (cross) and 
simulated reference temperature maps Tref(r) (solid line), la error bars are determined by 
applying the Monte-Carlo approach introduced in Schuecker et al. (2004). 
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Fig. 7. — Turbulence velocities Wturb as a function of turbulence length scales I, which are 
plotted for the radius of 100 (solid), 125 (dash), 150 (dash dot), 175 (dot), and 200 (thick 
solid) kpc, respectively. The vertical lines mark the turnover points obtained in the 
wavelet spectra (§4.1 and Fig. 5). 



